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ABSTRACT 

To first order in the deviation from scale invariance the inflationary potential 
and its first two derivatives can be expressed in terms of the spectral indices 
of the scalar and tensor perturbations, n and n T} and their contributions to 
the variance of the quadrupole CBR temperature anisotropy, S and T. In ad- 
dition, there is a “consistency relation” between these quantities: n T = -IT. 
We derive the second-order expressions for the inflationary potential and ? its 
first two derivatives and the first-order expression for its third derivative, in 
terms of n, n T , S, T, and dn/d In A. We also obtain the second-order con- 
sistency relation, n r = -if [1 + O.llf + 0.15(n - 1)J. As an example we 
consider the exponential potential, the only known case where exact analytic 
solutions for the perturbation spectra exist. We reconstruct the potential via 
Taylor expansion (with coefficients calculated at both first and second order), 
and introduce the Pade approximant as a greatly improved alternative. 
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1 Introduction 


In inflationary models quantum fluctuations excited on very small length 
scales (~ Hj ~ 10" 23 cm) are stretched to astrophysical scales (~ 10 25 cm) 
by the tremendous growth of the scale factor during inflation (Hi is the 
value of Hubble parameter during inflation) [1]. This results in almost scale- 
mvanant spectra of scalar (density) [2] and tensor (gravitational-wave) [3] 
metric perturbations. Together with the prediction of a spatially-flat Uni- 
verse they provide the means for testing the inflationary paradigm. The ten- 
sor fluctuations lead to cosmic background radiation (CBR) anisotropy and 
a stochastic background of gravitational waves with wavelengths from about 
1 km to over 10 4 Mpc. The scalar fluctuations also lead to CBR anisotropy 
and seed the formation of structure in the Universe. 

The amplitudes and spectral indices of the metric fluctuations can be ex- 
pressed in terms of the inflationary potential and its derivatives, evaluated 
at the value of the scalar field when astrophysically interesting scales crossed 
outside the horizon during inflation (from galactic scales to the presently ob- 
servable Universe, corresponding to the eight e-foldings about 50 e-folds or 
so before the end of inflation). Techniques have been developed for relating 
the scalar and tensor spectra to the potential and its derivatives in an ex- 
pansion whose small parameter is the deviation from scale invariance [4]. In 
particular, to first order in the deviation from scale invariance the spectral 
indices and the power spectra of the fluctuations today can be written as [5] 1 
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1 Several minor errors in Ref. [5] have been corrected here: the factors of Hi +n in 

Eqs. (A5, A7) should be 7f 0 4 ; the factor of l.l(n-l) in Eq. (A8) is more precisely 1.3(n- 1)- 

the factor of 1 .2ut in Eq, (A 14) is more precisely 1.4nj’. 
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Here k is the comoving wavenumber, x = m n V' /V measures the steepness 
of the potential, prime denotes derivative with respect to the scalar field that 
drives inflation, subscript 50 indicates that the quantity is to be evaluated 
50 e-folds before the end of inflation, 2 m Pl = 1.22 x 10 19 GeV is the Planck 
mass, H 0 is the present value of the Hubble constant, r 0 ~ 2Hq 1 is the 
present conformal age of the Universe, and 7 ~ 0.577 is Euler’s constant. 
Scale-invariant metric perturbations correspond to (n — 1) = = o. The 

functions T{k) and T T (k) are the transfer functions for scalar [7] and tensor [8] 
metric perturbations respectively; for kr 0 < 100, both T(k) and T T (k) 1. 

From these expressions the consequences of the scalar and tensor metric 
fluctuations may be computed. In particular, the contributions to the vari- 
ance of the angular power spectrum of the CBR anisotropy on large angular 

scales (/ « 200) that arise predominantly due to the Sachs- Wolfe effect are 
given by [8j 
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where r LSS ~ r 0 /( 1 + z LSS ) 1/2 ~ r 0 / 35 is the conformal age at last scattering 
(^lss — 1100) and ji is the spherical Bessel function of order 1. (We note 
in passing that both expressions are based upon the approximation that the 
Universe is matter-dominated at last-scattering; the small contribution of 
radiation, about 10%-20%, leads to corrections that are of the same order [8] 
and would have to be included in a more accurate treatment when the data 
so demand.) 

The point about which the potential is expanded is in principle arbitrary. However, 
the spectral indices n and n T can only plausibly be measured on scales from IMpc - 
10 Mpc and S and T depend upon perturbations on these same scales, so it makes sense 
to choose the expansion point to correspond to when these scales crossed outside the 
horizon during inflation; in addition, by taking A^to = 1 several expressions simplify. 
The precise number of e-folds before the end of inflation when these scales crossed outside 
the horizon depends logarithmically upon the energy scale of inflation and the reheat 
temperature, see Refs. [4, 6]; for the sake of definiteness we take this number to be 50, 
which can easily be changed to the correct value for a given inflationary model. 
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The contribution of scalar and tensor metric perturbations to the variance 
o the quadrupole CBR anisotropy can be computed numerically [5] 

5ee 5([qfJ 2 ) 

47T 
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47T 

where the dependence upon (n - 1) and n T is given to first-order. In eval- 
uatmg these expressions the effect of transfer functions is negligible as the 
integrals are dominated by kr a ~ 2. For simplicity, following footnote 2 
we henceforth omit factors of feoro) 1 -" and (* 50 ro)-^; they are easily re- 
inserted if needed. 
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1.1 First-order reconstruction 


Since 5, T, n, and n T are expressed in terms of the potential and its first two 
derivatives, one can invert the expressions to solve for the potential and its 
first two derivatives m terms of S', T, n, and n T plus a “consistency relation.” 
Those expressions are [9] 
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In the second expressions for the potential and its first two derivatives we 
have used the consistency relation to express n T in terms of Z as £ should 
be easier to measure [10], Note that the sign of V' cannot be determined as 
it can be changed by a field redefinition 0 - -<fi, though a specific choice 
here determines the signs of various later expressions. 
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This procedure actually generates the full second-order term for V 50 , while 
the other expressions are first-order. Our nomenclature and the overall re- 
construction strategy are explained and clarified below. 

1.2 The perturbative reconstruction strategy 

Before delving into the algebra of the more complicated second-order expres- 
sions, let us orient the reader with an overview of the perturbative recon- 
struction process [6, 9]. The goal of the process is to express the inflationary 
potential and its derivatives, 4° for l = 0, 1 , 2 , • • in an expansion whose 
terms are powers of observables, T, §, n, n T , dn/d In A, and so on. The 
choice of the set of observables is of course not unique. From these, the po- 
tential can be reconstructed over some finite interval via a series expansion. 
While one can hope to learn about the potential over the interval that affects 
astrophysical scales, it is probably not realistic to hope to learn much about 
the potential globally without some additional a priori knowledge (e.g., the 
functional form of the potential). 3 

The formal expansion parameter is the deviation from exact scale invari- 
ance, which can be expressed as (n - 1); the other observables we shall use 
are: n T ~ 0[(n - 1)], | ~ 0[(n - 1)], and dn/dlnX ~ C[(n - l) 2 ]. (In 
the scale-invariant limit, V{<j>) = const, and n = 1, n T = 0, | = 0, and 

dn/dlnX = 0.) The expression for 4° begins at order (n - l) l / 2 T, with 
higher-order terms (n - 1)*+'/’ T, k = 1, 2, 3, • • The expressions for higher 
derivatives of the potential are not only of higher order in (n - 1), but also 
involve more observables. For example, to lowest order the expression for 
V^o only involves T; that for V' 0 involves T and that for V& involves T, 
5 ’ and (n - 1), and, as we shall see, that for V" involves T, (n - 1), 
and dn/d In A. (We have not included n T in any of the lists, assuming im- 
plicitly that it is expressed in terms of the other observables by means of a 
consistency relation, as we now discuss.) Likewise, the expression for a given 
derivative involves additional observables as one goes to higher and higher 
order. 


The one possible exception involves the accurate measurement of the stochastic 
background of gravitational waves on scales from 1km to 3000 Mpc (corresponding to 
N ~ ® ~ 50) in which case the inflationary potential could be mapped out directly since 
the amplitude of the tensor perturbation on a given scale is related to the value of the 
potential. 
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As a matter of nomenclature, we shall refer to the lowest-order term for 
any or for the consistency relation as “first-order,” the next as “second- 
order,” and so on. The four expressions from Ref. [9] given above are respec- 
tively the second-order expression for R 5 o, the first-order expressions for V' 0 
and V” 0 , and the first-order consistency relation. In the next Section we shall 
use second-order results from Ref. [11] to derive the second-order expressions 
tor V50, Vq Q} and the consistency relation, as well as the first-order expression 
for R 5 o'. Both the second-order expression for and the first-order ex- 
pression for Vg'o involve the derivative of the scalar spectral index, dn/dlnX, 
which is likely to be very difiicult to measure. Higher-order expressions for 
the potential and its first three derivatives, as well as the fourth derivative, 
will involve yet another even less accessible observable, and thus there is little 
motivation at present for proceeding further in the perturbative expansion. 

An important feature of reconstruction is that the problem is overdeter- 
mined; specifically, a set of M > 3 observables can be expressed in terms 
of the potential and its first M — 2 derivatives. This implies a “consistency 
relation,” which, for increasing M, contains terms of higher and higher order. 
The lowest-order consistency equation, n T = -±f , has been much discussed 
(e.g., m Ref. [4]) and arises through Eqs. (1), (7) and (8) which express n T , 

and T in terms of V 50 and Vg 0 . 

Calculating higher derivatives alone, while keeping the calculation of each 
derivative to lowest order, does not lead to the correct second-order term in 
the consistency equation, and nor does calculating the second-order correc- 
tions to the derivatives present. One must systematically do both. The 
second-order version of the consistency equation is obtained by calculating 
the potential, its derivative and Eq. (1) to a higher order. Adding an extra 
order to the calculation of Vg 0 adds a new observable, (n — 1), which will 
appear in the consistency equation at second order. To account for there 
being still only a single consistency equation, there must be a new equation, 
and because (n- 1) has only entered at second-order in V' 0 , we only need the 
first-order equation for V^. The second-order consistency equation, which 
we calculate m this paper, therefore relates n r , f and (n - 1), with the last 
only appearing as a second-order correction. Were one to desire a calculation 
to yet higher order, the same pattern would persist; each existing derivative 
must be calculated to one extra order and the next derivative to lowest or- 
der, introducing a new observable. This will generate next-order terms in 
the consistency equation with the new observable appearing at that order. 
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However, this presently cannot be done as third-order expressions for V 50 and 
M>o have not been calculated. 

2 Second-order Reconstruction Reduced to 
Practice 

The reconstruction equations for the scalar potential and its first two deriva- 
tives, evaluated to second-order, are given in Ref. [ 11 ], though not in terms 
of cosmological observables. They are given in terms of the perturbation 
amplitudes A 2 G and A 2 S . Very roughly, As is the horizon-crossing amplitude 
of the density perturbation on a given scale and A G is the horizon-crossing 
amplitude of the tensor perturbation (in the Appendix we provide some re- 
lations between notation used in that paper and this one.) Our purpose here 
is to express these second-order expressions for the potential and its first two 
derivatives in terms of the measurable quantities n, dn/dhi\, n T , 5, and T. 

The amplitudes A 2 S and A% are related to the observables 5, T, n T and 
n by: 

A% = 0.70(1 - 1.3n r )T, A% = 9.6[1 - 1.15(n - 1 )]S, (13) 

where the (n - 1 ) and n T dependencies have been found by evaluating the 
Sachs- Wolfe integrals numerically. Both expressions are accurate to second- 
order. 

Before deriving second-order expressions for the potential and its deriva- 
tives, we calculate the second-order version of the consistency relation. It is 
obtained from Eq. (2.9) of Ref. [ 11 ], 

n T A% 

~T = [ 1 + 3e “ 2? ?]> (14) 

where to the required order the slow-roll parameters e and tj (defined in the 
Appendix) are given by 

e = ~n T / 2 , n = (n - l )/2 - n T . (15) 

This gives a simple and very useful relation for A G /A 2 S) 

A 2 

~£2 = ~0.5n T [1 — 0.5n T + 1.0(n — 1)] . (16) 
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Substituting into Eq. (13), we find the second-order consistency relation 
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- 0.8n r + 0.15(n - 1)J , 
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or 

T 

= ~7ux [l + — 0.15(n — 1)] . 

To the required order we can use the first-order truncation n T 
the brackets, thereby obtaining an alternative form, 


(18) 


— inside 



r T 

1 4- 0.11— + 0.15(n — 1) , 


(19) 


where n T is given in terms of the more accessible quantities (n - 1) and T. 

Independent measurements of n, n T and ~ provide a powerful test of the 
inflationary hypothesis; in the space of these parameters inflationary models 
must lie on the surface defined by Eq. (18). In Figure 1 we illustrate the 
inflationary surface both without and with second-order corrections. The 
second-order corrections break the degeneracy in the (n - 1) direction, as 
well as typically reducing J viewed as a function of n T and (n- 1). However, 
the portions of the surface that feature large corrections are not favored 
by present cosmological data, and further, are susceptible to higher-order 
corrections. (Indeed, well away from scale invariance the surface would be 
noticeably different even just using Eq. (19) instead of Eq. (18), which differ 
by third-order terms.) 

Obtaining the reconstruction equations is simply a matter of substituting 
into Eqs. (3.4), (3.6) and (3.15) of Ref. [11] for V, V' and V" respectively. 
We give two alternative forms for each, the first using n T and the second 
substituting | for n T using the second-order consistency equation. They are 
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+0 .6n T (n - 1) - 0.2 (n - l) 2 - 1.1 — T 

dlnAj ’ 

= 21 (n - 1) + 0.43^ + 0.073 (|) 

— 0.015^(ra - 1) - 0.2(n - l) 2 - l.l-±-| T. (22) 
° dlnAJ v ' 

These expressions are accurate to second-order. Naturally, they agree with 
the first-order expressions given earlier. 

Though no expression is given in Ref. [11] for V'", by using the lowest- 
order expressions for e, 77, and a third slow-roll parameter f, and Eq. (3.13) 
which relates the three to dn / din A, one can obtain the first-order expression, 


V 5o/ m P\ = ±104 yj-nr - _ g _j_ 4 / n _ y 

n T v ' ’ 


±104i/~ [?^ lnA 


where the overall sign is to be the same as that of V' , The second-order term 
would require yet another observable. As remarked in Ref. [11], even this 
first-order expression features the rate of change of the scalar spectral index, 
which is likely to be very difficult to measure. Realistically then, in the near 
term only the value of the potential and its first two derivatives are likely to 
be accessible to accurate determination. 

The final step in reconstructing the potential is to use d^/dN, N being 
the number of e-foldings, to find the range of <j> that corresponds to the eight 
or so e-foldings of inflation relevant for astrophysics. In effect, this introduces 
an additional small parameter, v = (50 - A r )/50 < 8/50, where TV ~ 42 is 
the number of e-folds before the end of inflation when the smallest scale 
of interest crossed outside the horizon. The smallness of this parameter is 
essential if the approximation of the spectra by power-laws is to be valid. 

To proceed, we may simply carry out a Taylor expansion of <j> about 0 5O , 
to whatever order we believe is appropriate, 


(f>N ~ 050 = (N — 50) 


— + ~(N - 50) 2 --- 

dN *<, 2( } 
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We can easily substitute for these derivatives, 
exact formula 




47T 


using as a starting point the 

(25) 


which, along with dN/dt — - H , yields the relation from which the Taylor 
coefficients may be calculated, 


d<f> _ m pj H' 

dN = "i7 ( 26 ) 

To get a given coefficient in the Taylor expansion for <j > N , one simply 
calculates d l cf)/dN l expanding to the desired order in the deviation from scale 
invariance. 4 For example, taking only the first term in the (p N expansion and 
working to first-order one obtains 

(f)N - 050 ^ ±^~y/-n T (N — 50), (27) 

where again the overall sign is the same as that of V\ Next, we give the first 
coefficient in the cj) N expansion to second-order and the second coefficient in 
the 07v expansion to first-order only, 

0JV-05O = ± -~y/-n T [1 + 0.1 n T + 0.1(n - 1)] (N - 50) 

Tfl pj 

_ _ nr) (iV - 5 °) 2 + ’ ' ' - ( 28 ) 

with both signs again agreeing with that of V' . 


3 Techniques and an Example 

Before going on to specifics, let us again consider perturbative reconstruction 
from the larger perspective. In constructing the Taylor series for the value of 

4 Note this procedure differs slightly from that in Ref.- [9], where d<f>/dN was expanded 
linearly about 0 5O and (j> N was solved for exactly, cf. Eq. (8). This results in an exponential 
whose expansion picks up the (N- 50) and (iV-50) 2 terms correctly to lowest order in the 
deviation from scale invariance, though not the higher-order terms in the (N - 50) term 
which would require higher-order terms in the expansion of d<f>/dN. There is an overall 
sign error in Eq. (8) of Ref. [9]. 
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the potential at a given point the l-th term is given by V$((f) - (t> 50 ) l /l\. As 
mentioned earlier starts at order (n - 1 ) l / 2 T; from Eq. (28) we see that 
(* - *»)' -der (n- - 1) //2 AA^, where AN = A - 50 is around eight for 
the region of the potential corresponding to astrophysically accessible scales. 
Thus, the l-th term in the Taylor expansion to lowest order is 

V 5o(<t> ~ foo) l /H ~ 0[(n - l) 1 AN 1 T] 4 . (29) 

In order that the series be clearly convergent, the deviation from scale in- 
variant (n - 1) should be less than about A TV -1 ~ 0.1. In the case that 
(n - 1)AN ~ 0(1), as in our example below, Eq. (29) suggests that the 
expression for each derivative be expanded to the same order in the devia- 
tion from scale invariance. On the other hand, if (n - 1) < 0.1, then clearly 
the higher-derivative terms in the Taylor expansion are less important, and 
a sound case could be made for expanding the lower derivatives of the po- 
tential to higher order in (n - 1). However, because of the existence of two 
expansion parameters, AN and (n - 1), that are a priori unrelated, in the 
general case there is no clear prescription. In our example below, where 
(n - l)AiV ~ 0(1), we will explore several possibilities. 

3.1 Expansion techniques 

Given the value of the potential and its first two or three derivatives at a 
point and the <f) N relation just obtained, one can reconstruct the potential 
on the observation ally relevant scales (i.e., N ~ 42 - 50). The standard 
technique used previously is the Taylor expansion 

= ^50 + I 50 O - 05o) + - </> 5 o ) 2 4 (30) 

For many situations this is perfectly fine (e.g., when tit and n — 1 are small, 
see Ref. [9]). However, if the range of eight or so e-foldings corresponds to a 
large range in 0 the convergence may not be very good because of the abrupt 
truncation of the Taylor series. Specifically, for large (0 - </> 50 ) the shape of 
the reconstructed potential is dictated, rightly or wrongly, by the last term 
in the expansion (quadratic or cubic). 

An often used alternative is the Pade approximant [12], which can be 
generated directly from a truncated power series. For a power series that 
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extends to order N, the Pade approximants are quotients of two polynomi- 
als of order L (numerator) and M (denominator) denoted by [L,M], where 
L + M — N. By construction, the expansion of [L, M] matches that of the 
power series to order N , but of course is not truncated. Very often, Pade ap- 
proximants provide a very good approximation over a wider range of values 
than the Taylor series from which they are derived; they in some way encode 
better estimates of the higher-order terms than does truncation. If we trun- 
cate the Taylor series at the second derivative, then the associated diagonal 
Pade approximant [1, 1] is a ratio of two first-order polynomials given by 5 


m = 


do a i cf) 

1 + bi<f) 


(31) 


with 

a 0 = V 50 ; h = -V!' 0 /2V' 0 ; a, = V' 0 - V 5Q V'' 0 / 2V' 0 . (32) 

As we shall now illustrate by specific example, Pade approximants have a lot 
to offer when the Taylor series proves a poor approximation. 


3.2 Reconstructing an exponential potential 

A useful testing ground for reconstruction is the exponential potential, the 
only known case where the perturbation spectra can be derived exactly an- 
alytically [13]. For the potential 


l ' w = F,exp (-fi)' ( 33 ) 

the scale factor grows exactly as t p . Compared with the lowest order expres- 
sions, the amplitudes A and A T , or A% and A%, are both multiplied by the 
same p-dependent factor R 2 (p ), where 

R(p) = 2 V(r-D (1 _ 1/p) P/(P- 1 ) _ (34) 

The [2, 0] approximant is just the truncated Taylor series; in addition to simplicity, 
there is some motivation for using the diagonal approximant rather than the [0, 2] ap- 
proximant as it is asymptotically constant, consistent with the flatness of inflationary 
potentials. 
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where T(- • •) is the usual gamma function. Both scalar and tensor spectra 
are exact power laws with spectral indices (n — 1) = ut ~ — 2/(p — 1). The 
scalar-field solution is characterized by 

• / p rupj d(b mpi 

^ = V^~ ; dN = ~^ ^) = ^ rue exp(2(lV-50)/ P ]. (35) 

The expressions for T and S can be obtained exactly by integrating Eqs. (5,6). 


5 = 2.2/(n)fl 2 (p)^»_ 

m pi x 5o 
T/true 

T = 0.61 g(n T )R 2 (p)\-^ 

mh 


(36) 

(37) 


where the numerical factors f(n) = 1 + 1.15(n - 1) + • • • and g(n T ) = 1 + 

l.3n T + ••• arise from the n, n T dependence of the Sachs- Wolfe integrals cf 
Eqs. (5, 6). 

We are now ready to carry out an array of reconstruction methods. Be- 
cause we are using exact expressions to generate the spectra, this procedure 
is more ambitious, and more realistic, than those attempted thus far [6, 9], 
where the trial spectra were produced using the slow-roll approximation. For 
the general inflationary potential, exact results are not known, and so this 
procedure is not possible. However, our method here should give a more 
realistic estimate of inherent errors even in the general case. 

There are two distinct types of error. The first is error in the value of the 
potential at </> 50 , due to third-order and higher terms. By substituting the 

expression for T in Eq. (37) into Eq. (9) or (20) for V 50 we can compute that 
error: 


- g{n T )( 1 - 1 An T )R(p) 2 . ( 38 ) 

The second error involves the shape of the potential, which depends on the 
ability of the chosen expansion to match the potential over the eight inter- 
esting e-foldings. 

We have chosen as a specific, example an exponential potential with p = 
43/3. We did so because this leads to about the largest departure from scale 
invariance that can still be regarded as observationally viable, (n — 1) = 
n T = -0.15 and | ~ 1, and thus realistically represents the most challenging 
example of reconstruction. The exact potential is shown in Fig. 2 along with 
the results of five different reconstructions. 
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To begin, consider the error in estimating Vj 0 rue ; g(n T = -0.15) = 0.824 
and so V 5 o/Vf/ ue ~ 0.95, a modest 5% due to the neglected higher-order 
terms. As we always include the second-order term in R 50 , the error is the 
same in every method we look at. Had the first-order expression for R 50 been 
used instead, corresponding to the neglect of the factor of (1 - l.4n r ) i n 
Eq. (9), then the underestimation would have been about 20%. 

Let us now consider the shape, which we note depends on T and S only 
through their ratio. The important distinction between different methods is 
the difference in required input data; methods needing only n and £ have 
the advantage of depending only on the information that is easiest to obtain. 
Requiring dn/d In A in addition, while offering more accuracy, is setting a 
much trickier observational task. 

As a starting point, let us take the equations derived in Ref. [9], which 
are primarily first-order though they include the second-order correction to 
V 50 , cf. Eqs. (9-11). In this extreme example, the quadratic Taylor series 
based upon this does a bad job of approximating the shape of the potential, 
as it turns upward for large {<j> - (f> 50 ) due to the truncation at the (<b — d ^) 2 
term (see Fig. 2). 

If we now require knowledge of dn/dln A, the Taylor series approach can 
be improved in two ways. We can now take R 50 , R 5 ' 0 , and R" to second- 
order; however, the improvement is rather minimal. Alternatively, we can 
stick to first-order expressions, but include the V 5 'g cubic term. Again the 
improvement is modest, though at least the unwanted minimum has been 
eliminated. One could go further and take R so , V' 0 , and R" to second- 
order and V 5 q to lowest order, which we haven’t illustrated, again seeing 
only modest gains for the increased observational requirement. 

The Taylor series having been unimpressive, let us progress in a different 
direction. With only n and £, as an alternative to the Taylor series one can 
construct the Pade approximant based upon it, taking R 50 to second-order 
and R 5 o and V" Q to first-order. This represents a substantial gain on the Tay- 
lor series to that order without requiring any additional input information. 
With this minimal information, it is a much better method. Reintroducing 
dn/d In A allows this method to be extended to second-order, where the re- 
production of the shape of the potential is excellent. To include the third 
derivative term would necessitate a more complicated (non-diagonal) Pade 
approximant, which doesn’t seem warranted at the moment. 

What is the upshot of this comparison? Recalling that we have chosen 
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an example with extreme deviation from scale-invariance, the second-order 
corrections are reassuringly small and only improve the shape of the recon- 
structed potential slightly. The addition of the third derivative term in the 
Taylor series gives a slightly more significant improvement, but at the price 
of its dependence upon dn/d In A even at lowest order. The most remarkable 
improvement involves the use of Pade approximants. Even without knowl- 
edge of dn/d In A the shape of the potential is reproduced far better than 
with the higher-order Taylor series which does require that knowledge. As 
noted previously, the improvement results from the fact that the Pade ap- 
proximant is not truncated; further, even in situations where truncation of 
the Taylor series does not lead to problems, the Pade approximant still proves 
valuable as its Taylor expansion coincides with that of the original expan- 
sion. We therefore conclude that Pade approximants provide a significant 
improvement in the perturbative reconstruction of the inflationary potential. 


4 Discussion 

By presenting the second-order reconstruction equations directly in terms 
of observables, we have been able to assemble and to compare an array of 
different perturbative reconstruction techniques based upon cosmological ob- 
servables. Our work extends previous work in several important ways. 

The most interesting result is the introduction of the Pade approximant 
as an alternative to the Taylor series in perturbative reconstruction. It can 
be obtained from a Taylor series regardless of the order (in the deviation 
from scale invariance) to which the coefficients of the Taylor series has been 
obtained. In our worked example, the improvement in reproducing the shape 
of the potential as compared to the Taylor series is striking, especially con- 
sidering that no extra observables are required. 

We have shown that the second-order corrections to the Taylor series 
coefficients are generally small, and that those for F 50 and V 5 ' 0 only depend 
upon the same quantities as the first-order expressions (5, T, and n). The 
corrections to V/' 0 however require a new observable such as dn/d In A, and 
by deriving for the first time an explicit expression we have confirmed that 
even the lowest-order term in V/// requires this challenging observable. 

Finally, one of the most important aspects of reconstruction is that it 
is over determined: Any set of cosmological observables supplies degenerate 
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information regarding the potential and its derivatives, thereby providing 
an important consistency check. In particular, the tensor spectral index 
can be expressed to second-order in terms of S, T, and n by the relation: 

Ut = ~7sl 1 + °- n | + 0.15(n - 1)]. In cases that are observationally viable, 
the second-order corrections are small. 
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Appendix: Some Relations between Notation 


For the convenience of the reader, we summarise here some relations between 
the notation used here and that in Ref. [11], from which several important 
results were taken. In that paper, the spectra A s and A G were defined so 
as to include any scale-dependence within them, i.e., they are functions of k. 
In the perturbative reconstruction regime, where for most results the spectra 
can be approximated by power laws, these are related at lowest-order to the 
amplitudes A and A T in this paper, which are just numbers, by 

?7T 2 

Mk/ks or 1 = ~Al(k), (39) 

A T (k/k 50 ) nT = 2 A%(k). (40) 

In making the connection, note that 87 t ~ 25. For the higher-order terms, 
the scalar field kinetic energy must be accounted for in translating between 
H and V, which means these relations break down at higher order. 

In Ref. [11], slow-roll parameters e and r) are introduced, 


e = H!k(K ) 2 n = 

4 TT \H ) ' V 47 T H' 


(41) 
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which are again in general ^-dependent. As indicated in Section 2 of the 
present paper, they can be related to the spectral indices to various orders, 
e and r] being of the same order in perturbation theory as (n — 1) and n T . 
To lowest-order they are constant, corresponding to power-law spectra. At 

lowest-order e = 167nr 2 , but once more higher order corrections break this 
relation. 
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Figure Captions 


Figure 1: The consistency plane for inflation in n — n T — | space, the flat 
surface being the lowest-order result and the curved one incorporating the 
second-order corrections, given by Eq. (18). 

Figure 2: An array of different reconstructions of an exponential potential 
with (n — 1) — n T = -0.15 (p = 43/3). The longer dotted line indicates 
the exact potential. The three different line styles correspond to three differ- 
ent reconstruction strategies; solid is Taylor series truncated at — 0 5o ) 2 , 
dashed is Taylor series truncated at - <^ 0 ) 3 an d dash-dotted is the Padi 
approximant based on the former of these. The upper line of a given style 
uses coefficients to first-order in the deviation from scale invariance (save V 50 , 
which is always second-order), while the lower, where plotted, is second-order 
in all coefficients. The length of the curves corresponds to eight e-foldings. 
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